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CN ! Abstract 

\ We present the calculation of the vh'tual- and bremsstrahlung corrections of 0{as) to the matrix 

• elements {dl'^ l!.^\Oi\h) . This is the missing piece in the full next-to-next-to- leading logarithmic 

(NNLL) results for various observables associated with the process B Xf^i~^l~ , like the branching 
J> ' ratio, the CP-rate asymmetry and the forward-backward asymmetry. This paper is an extension 

m 

^3 ' of analogous calculations done by some of us for the process B Xsi^f~ . As the contributions 

■ of the diagrams induced by the four-quark operators O" and with a u-quark running in the 



o 



quark loop are strongly CKM suppressed, they were omitted in the analysis B ^ Xgi^i . This 



' is no longer possible for B X^l^i , as the corresponding contributions are not suppressed. The 

Q^' main new work therefore consists of calculating the 0{as) corrections to {d l~ \0i 2\b) ■ In this 

paper we restrict ourselves to the range U.U5 < s/m^ < U.25 (s is the invariant mass of the lepton 
pair), which lies above the p- and w-resonances and below the JZ-i/^-resonance. We present the 

X" 

' analytic results for the mentioned observables related to the process B — > X^vi as expansions 



in the small parameters s = s/rn^, z = ml/rn^ and s/(4m^). In the phenomenological analysis at 
the end of the paper we discuss the impact of the NNLL corrections on the observables mentioned 
above. 
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I. INTRODUCTION 



It is well-known that various observables associated with inclusive rare S-decays like 
B — > Xs^dl s-nd B — > sensitively depend on potential new physics contributions. 

But even in the absence of new physics these observables are important, because they provide 
checks on the one-loop structure of the Standard Model (SM) theory and can be used to 
gain information on the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements Vts and Vtd, 
which are difficult to measure directly. 

At present, a lot of data already exists on BR{B — > Xg^) [1-7] and on BR{B — > Xsi~^i~) 
[8-10] and it is expected that in the future also data on the CKM suppressed counterparts, 
i.e. on BR(S — > Xa'f) and on BR{B Xai'^£~) will become available. The same holds for 
experimental information on additional observables, like CP-rate asymmetries or forward- 
backward asymmetries in the decays B — > Xsj£^£~ . 

In order to fully exploit and interpret the experimental data, it is obvious that precise 
calculations in the SM (or certain extensions thereof) are needed. The main problem in the 
theoretical description of the decay B — > Xsi'^£~ is due to the long-distance contributions 
induced by cc resonant states and in principle also by uu resonant states. The latter are, 
however, strongly CKM suppressed. This suppression is not present in the case oi B ^ 
XJ.^£r ^ as the CKM factors involved in the contributions from cc and uu resonant states 
are of the same order. When the invariant mass y/s of the lepton pair is close to the 
mass of a resonance, only model dependent predictions for such long distance contributions 
are available at present. It is therefore unclear whether the theoretical uncertainty can be 
reduced to less than ±20% when integrating over these domains [11]. 

However, restricting ^/s to a region below the cc resonances, the long distance effects in 
B Xgi^i^ are under control. The same is true for B X^i^i^ when choosing a region 
of y/s which is below the J/ip- and above the p, o^-resonance regions. It turns out that in 
those ranges of ^/s the corrections to the pure perturbative picture can be analyzed within 
the heavy quark effective theory (HQET). In particular, all available studies indicate that 
for the region 0.05 < s — s/ml < 0.25 the non-perturbative effects are below 10% [12-17]. 
Consequently, observables like differential decay rates, forward-backward asymmetries and 
CP-rate asymmetries for B Xs^^i'^i^ can be precisely predicted in this region of ^/s using 
renormalization group improved perturbation theory. It was pointed out in the literature 
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that the differential decay rate and the forward-backward asymmetry in S — > Xsi'^i~ are 
particularly sensitive to new physics in this kinematical window [18-20]. 

In the context of the SM there exist computations of next-to-leading logarithmic (NLL) 
QCD corrections to the branching ratios for B — > Xg^ [21-28] and B X^'j and the 
corresponding CP-rate asymmetries [29-31]. Next-to-next-to-leading logarithmic (NNLL) 
QCD corrections to the branching ratio [32-36] and the forward-backward asymmetry in 
B — > Xsi'^£~ are also available [37-40]. For a recent review see e.g. [41]. 

The corresponding NNLL results for the process B — > Xdi'^i~ are missing, however. The 
aim of the present paper is to close this gap. The main difference between the calculations 
for B Xs£~^£~ and B X^i^i^ lies in the contributions of the current-current operators. 
In the existing NNLL calculations of 5 — > X^i'^i" only those associated with and O2 
were included at the two-loop level because those induced by and O2 are strongly CKM 
suppressed (see Section II for the definition of the operators 0^'2). For B — > Xdi'^i~ the 
contributions generated by O" and O2 are no longer CKM suppressed and have to be taken 
into account as well. At first sight, it seems that the two-loop matrix elements of O" ^^<i 
O2 can be straightforwardly obtained from those of and O2 by simply taking the limit 
rric — >■ 0. This is, however, not possible for some of the diagrams in Fig. 1, because the two- 
loop matrix elements of 0^ and O2 were derived by doing various expansions. In particular, 
one of the expansion parameters is s/(4m^), which is formally <^ 1 when restricting y/s to the 
window discussed above. Obviously, the analogous quantity for the li-quark contribution, 
s/(4m^), cannot be used as an expansion parameter, which implies that genuinely new 
calculations for the w-quark contributions are needed. As discussed in Section III, the 
calculations of certain diagrams associated with even more involved than those 

associated with 0^2- To derive the new results, we used dimension-shifting techniques in 
order to reduce certain tensor integrals to scalar ones and integration-by-parts techniques 
to further simplify the scalar integrals [42, 43]. 

As the main emphasis of this paper is the derivation of the matrix elements {d£'^i~\Oi 2\b) 
at order as, we keep the phenomenological analysis relatively short. In particular, we do 
not take into account power corrections, but merely illustrate how the NNLL contributions 
modify the scale dependences of the branching ratio, the forward-backward asymmetry and 
the CP-rate asymmetry. A more detailed phenomenology, including power corrections, will 
be presented elsewhere. 
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The paper is organized as follows: In Section II we present the effective Hamiltonian for 
the decay b — > X^l^l' . Section III is devoted to the virtual 0{as) corrections to the opera- 
tors O"'*^ and O"''^. Subsequently, Section IV presents the corresponding contributions to the 
form factors of the operators O7, Og, O9 and Oio- With these results at hand, we discuss in 
Section V the corrections to the decay width of -B — > Xd^^^~ . In Section VI we show some 
applications of our results. A summary of the paper is presented in Section VII. The appen- 
dices contain technical details about the performed calculation: Appendix A explains the 
dimension-shifting and integration-by-parts techniques. These techniques are then applied 
to the calculation of diagrams Id), which is presented in Appendix B. Appendix C outlines 
a procedure on how to calculate the evolution matrix for the Wilson coefficients as a power 
series in a^. Appendix D contains one-loop matrix elements needed in the calculation of 
the counterterms. Finally, in Appendix E we present the results for those bremsstrahlung 
contributions which are free of infrared and coUinear divergences. 

II. EFFECTIVE HAMILTONIAN 

The appropriate framework for studying QCD corrections to rare 5-decays in a systematic 
way is the effective Hamiltonian technique. For the specific decay channels B X^i^i^ 
and B Xiii~^i~ {£ = fi, e), the effective Hamiltonian is derived by integrating out the 
i-quark, the H^-boson and the Z°-boson. In the process B Xs£~^£~ , the appearing CKM 
combinations are A„, Ac and At, where A^ = VibV*^. Since |Au| is much smaller than |Ac| and 
I At I, it is a safe approximation to set A^ equal to zero. Using then the unitarity properties 
of the CKM matrix, the CKM dependence of the Hamiltonian can be written as a global 
factor \t. In the case of B Xfii~^i~, all three quantities C,i = yibVid = u,c,t) arc of the 
same order of magnitude. Therefore, as no approximation is possible, the CKM dependence 
does not globally factorize. The effective Hamiltonian reads 
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We choose the operator basis according to [32] : 

O3 = ((iL7M^L)E,(^7'^?), O4 = {dL'y^.T%L)Egi<l'^''T'^'i)^ 

(2) 

O7 = ^mb{dLa^''bR)F^,, Os = j-mb{dLa^''T%R)G';^,, 

O9 = |(47m^'l)E^(^V^), Oio = J(rf"L7A)E^(^r75^), 

where the subscripts L and R refer to left- and right-handed components of the fermion 
fields, respectively. 

The factors l/g^ in the definition of the operators O7, Og and Oio as well as the factor 
I/qs present in Og have been chosen by Misiak [44] in order to simplify the organization of 
the calculation. With these definitions, the one-loop anomalous dimensions [needed for a 
leading logarithmic (LL) calculation] of the operators are all proportional to g^, while 
two-loop anomalous dimensions [needed for a next-to-leading logarithmic (NLL) calculation] 
are proportional to g^, etc. 

After this important remark we now outline the principal steps which lead to a LL, NLL, 
and a NNLL prediction for the decay amphtude for 6 — > d i'^i" : 

1. A matching calculation between the full SM theory and the effective theory has to be 
performed in order to determine the Wilson coefficients Q at the high scale fiw ~ 
m.Wj'iTit- At this scale, the coefficients can be worked out in fixed order perturbation 
theory, i.e. they can be expanded in g^: 

Qinw) = Ci'\nw) + ^2Cl'\^^w) + ^^.CPi^iw) + 0{g!) . (3) 

At LL order, only Cf*^ are needed, at NLL order also Cf^\ etc. The coefficient Cy^^ 
was worked out in Refs. [23-25], while Cg and C)q were calculated in Ref. [32]. 

2. The renormalization group equation (RGE) has to be solved in order to get the Wilson 
coefficients at the low scale fib ~ rub. For this RGE step the anomalous dimension 

matrix 7(q;s), which can be expanded as 

.(0) ^ , ^(1) ^ , ^(2) / ^ ^ 



is required up to the term proportional to 7^^) when aiming at a NNLL calculation. 
After the matching step and the RGE evolution, the Wilson coefficients Cj can be 
decomposed into a LL, NLL and NNLL part according to 

C>M = Cr M + + + 0(<,») . (5) 

We stress at this point that the entries in 7^^^ which describe the three-loop mixings 
of the four-quark operators Oi — Oq into the operator O9 have been calculated only 
recently [33]. In order to include the impact of these new ingredients on the Wilson 
coefficient Cg^^h), wc had to reanalyze the RGE step. In Appendix C, we derive 
a practical formula for the evolution matrix U{iJ,b, ^w) at NNLL order, generalizing 
existing formulas at NLL order (see e.g. [45]). 

3. In order to get the decay amplitude, the matrix elements ((i£+£^|0j(/i;,)|6) have to be 
calculated. At LL precision, only the operator O9 contributes, as this operator is the 
only one which at the same time has a Wilson coefficient starting at lowest order and 
an explicit 1/gl factor in the definition. Hence, at NLL precision, QCD corrections 
(virtual and bremsstrahlung) to the matrix element of O9 are needed. They have 
been calculated in Refs. [44, 46]. At NLL precision, also the other operators start 
contributing, viz. Oj{fib) and Oio(a*6) contribute at trcc-lcvel and the four-quark 
operators Oi,...,Oq at one-loop level. Accordingly, QCD corrections to the latter 
matrix elements are needed for a NNLL prediction of the decay amplitude. 

The formally leading term ~ [1 / g1)Cf\iJ,}j) to the amplitude for h — > d£~^£~ is smaller 
than the NLL term ~ {l/gl)[g^J{167r^)]C^^\iXb) [47]. We adapt our systematics to the 
numerical situation and treat the sum of these two terms as a NLL contribution. This is, 
admittedly some abuse of language, because the decay amplitude then starts out with a 
term which is called NLL. 

As pointed out in step 3), 0{as) QCD corrections to the matrix elements {di~^i~\Oi{^b)\b) 
have to be calculated in order to obtain the NNLL prediction for the decay amplitude. In 
the present paper we systematically evaluate virtual corrections of order as to the matrix 
elements of Oi, O2, O7, Os, Og and Oiq. As the Wilson coefficients of the gluonic penguin 
operators O3, Oe are much smaller than those of Oi and O2, we neglect QCD corrections 
to their matrix elements. We also systematically include gluon bremsstrahlung corrections 
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b 0^,2 d b 0^,-2 d b 0^,^ d 




b) d) f) 

FIG. 1: Complete list of two- loop Feynman diagrams for b ^ d^* associated with the operators 
O"'*^ and Oa''^. The fermions (6-, d-, u- and c-quarks) are represented by solid lines, whereas the 
curly lines represent gluons. The circle-crosses denote the possible locations where the virtual 
photon (which then splits into a lepton pair) is emitted. 

to the matrix elements of the operators just mentioned. Some of these contributions contain 
infrared and coUinear singularities, which are canceled when combined with the virtual 
corrections. 

III. VIRTUAL 0{ots) CORRECTIONS TO THE MATRIX ELEMENTS 

{d£+e-\oT^^\b) 

In this section we present the calculation of the virtual 0{as) corrections to the matrix 
elements of the current-current operators O"''^ and 02 "^. Using the naive dimensional regular- 
ization scheme (NDR) in c? = 4 — 2 e dimensions, both ultraviolet and infrared singularities 
show up as 1/e" poles {n = 1,2). The ultraviolet singularities cancel after including the 
counterterms. CoUinear singularities are regularized by retaining a finite down quark mass 
rud- They are canceled together with the infrared singularities at the level of the decay width 
when taking the bremsstrahlung process b dt^l~g into account. We use the MS renormal- 
ization scheme, i.e. we introduce the renormalization scale in the form /Z^ = exp(7£;) / (4 tt) 
followed by minimal subtraction. The precise definition of the evanescent operators, which 
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is necessary to fully specify the renormalization scheme, will be given later. 

Gauge invariance implies that the QCD corrected matrix elements of the operators 
can be written as 

{den-\Ol\b) = ^;',^09)tree + {07)t.ee (i = 1, 2; g = U, c) , (6) 

where (09)tree a-nd (07)tree are the tree-level matrix elements of Og and O7, respectively. 
Equivalently, we may write 

{den-\0',\b) = [i^^:^(09)tree + F^l\dr),.^] , (7) 

where the operators O7 and Og are defined as 

07=^07, 09 = ^09. (8) 

We present the final results for the QCD corrected matrix elements in the form of Eq. (7). 
The full set of the diagrams contributing at 0{as) to the matrix elements 

Mf^{d£n-\0!\b) (9) 

is shown in Fig. 1. As indicated, the diagrams associated with O"''^ and 02'"^ are topologically 
identical. They differ only in the color structure. While the matrix elements of the operator 
Og"^ all involve the color structure 



there are two possible color structures for the corresponding diagrams of O^''^, viz 

Ti = ^ T^'T^T^'T^ and T2 = ^ jarpbrpbrpa 

a,b a,b 

The structure Ti appears in diagrams la)-d), and T2 enters diagrams le) and If). Using the 
relation 



we find that ri = C-ril and T2 = C^-jl, with 



= 47V^ = 4^2 



8 



Inserting A^^, = 3, the color factors are = |, = — | and C^^ = ^. The contributions 
from O"'"^ are obtained by multiplying those from Og''^ by the appropriate factors, i.e. by 
Cti/Cf = — I and Ct^/Cf = |, respectively. As the renormalized 0{as) contributions of 
the operators 01 and O2 are discussed in detail in Ref. [34], we only discuss the calculations 
of the contributions from and 0| to the individual form factors. 

The rest of this section is organized as follows: We discuss the calculations of the di- 
agrams la)-e) for the operators 0^2- Notice that all results are given as an expansion in 
the small quantity s — s/ml, where s is the invariant mass squared of the lepton pair, 
and that we keep only terms up to 0{s^). After deriving the counterterms that cancel the 
divergences of the diagrams mentioned above, we present the renormalized contributions to 
the form factors. We postpone the discussion of diagrams If) as it turns out to be more 
convenient to take them into account when discussing the virtual corrections to Og. 



A. Diagrams la) and b) 

The calculation of the contributions to F^"^^ and F^^^ from the diagrams in Figs, la) and 
lb) opposes no difficulties, as it can be performed by using the Mellin-Barnes approach [48] . 
Alternatively, one may get the results directly from the corresponding form factors of the 
b si~^i~ transition by taking the limit rric 0. The form factors associated with the 
diagrams in Fig. la) are given by 
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where 



Ls = ln(s) and Ly = ln( — 
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For the sum of the diagrams in Fig. lb) we find 
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B. Diagrams Id) 



The computation of the diagrams in Fig. Id) is by far the most comphcated piece in 
our entire calculation of the 0{as) corrections to the matrix element for b d£'^£~. After 
various unsuccessful attempts, we managed to obtain the result by using the dimension- 
shifting method [42] (see Appendix Al), combined with the method of partial integration 
(see Appendix A 2) . Since we want to include the details of the actual calculation, we relegate 
them to Appendix B. Here, we merely present the final results, viz the contributions to the 
form factors, which read 
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C. Diagrams Ic) 



The calculation of this diagram can be done in a very simple and efficient way. We add 
the two subdiagrams and integrate out the loop momentum of the virtual gluon. Next we 
integrate over the remaining loop momentum, being left with a four dimensional Feynman 
parameter integral. After introducing a single Mellin-Barnes representation of the occurring 
denominator, the parameter integrals can all be performed. At this level, the result contains 
Euler Beta-functions involving the Mellin-Barnes parameter. Finally, the Mellin-Barnes 
integral can be resolved applying the residue theorem, which naturally leads to an expansion 
in the parameter s. The contribution of the diagrams in Fig. Ic) to the form factors reads 
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^6 3 3 9 y ' ■ 3 
We also performed the calculation of this diagram in two different, more complicated ways, 
namely by 

• using the building block Jq,^ given in [34] and then introducing a double Mellin-Barnes 
representation, 

• using the dimension-shifting and integration-by-parts techniques as explained when 
discussing the computation of the diagrams in Fig. Id) (see also Appendices A 1 and 
A2). 

We found that all three calculations yield the same result and thus serve as an excellent 
check for the dimension-shifting approach and for the very complicated double Mellin-Barnes 
calculation. 
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D. Diagrams le) 



The diagrams in Fig. le) may again be solved in two ways. The first way is to use the 
large external momentum expansion technique [48]. The second possibility is to apply the 
dimension-shifting and integration-by-parts procedure [42] also for this diagram. We do 
without presenting the calculation and merely give the results for the contributions to the 
form factors. 
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E. 0{a.s) counterterms to {dl+l-\0'^l^\h) 

So far, we have calculated the two-loop matrix elements (d | Q 0|^| 6) (i — 1,2, 
q — u,c). As the operators mix under renormalization, there are additional contributions 
proportional to Cj. These counterterms arise from the matrix elements of the operators 

2 10 12 

J^^Z^.^OJ + Op+J^^^O. + J^^Z^.^OJ + Op, i = l,2, (17) 

j=l 3=2, j=\l 

where the operators Oi-Oio are given in Eq. (2). O^^" and O^^ evanescent operators, i.e. 
operators which vanish in c? = 4 dimensions. In principle, there is some freedom in the choice 
of the evanescent operators. However, as we want to combine our matrix elements with the 
Wilson coefficients calculated by Bobeth et al. [32], we have to use the same definitions: 

On = [dLl^lulaT^Ui) (uLYYrT%L) - 16 

0^2 = (^TmT.T.^^l) {uLrrrbL) - 16 , (18) 
OJi = (47M7.7a7^"ci) (cLrrrT^L) - 16 0J , 



The operator renormalization constants Zij = 6ij + 6Zij are of the form 

f 02 ^ 

Any ' e^*V ' (I^V'^'^e^*^ ' e 
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The coefficients a-^ needed for our calculation we take from Refs. [32, 34] and list them for 
i = 1,2 and j = 1, 12: 
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6 00 I 0000 -| 010 



(20) 



„12 _ __58_ „12 _ _64_ 22 _ 1168 
"17 243 ' 19 ~ 729 ' 19 ~ 243 ' 

(21) 

_12 _ 116 „12 _ 776 „22 _ 148 
"27 ~ 81 ' 29 ~ 243 ' "29 ~ 81 ' 

We denote the counterterm contributions to b ^ d t^t~ which are due to the mixing of 0\ 
or O2 into four-quark operators by F^^S-lAquax^i ^^"^ -^i'^u-l4quark- They can be extracted from 
the equation 

(22) 

where j runs over the four-quark operators. The operators are understood to be identified 
with Oj for j — 3, 4, 5, 6. As certain entries of a}^ are zero, only the one-loop matrix elements 
of O"''^, 02''^, O^'^, and are needed. In order to keep the presentation transparent, 
we relegate their explicit form to Appendix D. We do not repeat the renormalization of the 
01 and O2 contributions at this place and refer to [34]. 

There is a counterterm related to the two-loop mixing of O" (i = 1, 2) into O7, followed by 
taking the tree-level matrix element {d i~ {Ojlb) . Denoting the corresponding contribution 
to the counterterm form factors by F^^^J}^ and F^^^^^, we obtain 

FS2r--f, i^S7 = 0. (23) 

The counterterms which are related to the mixing of {i = 1,2) into O9 can be spht 
into two classes: The first class consists of the one-loop mixing O9, followed by taking 

the one-loop corrected matrix element of O9. It is obvious that this class contributes to 
the renormalization of diagram If), which we take into account when discussing the virtual 
corrections to O9. We proceed in the same way with the corresponding counterterm. 

The second class of counterterm contributions due to Of — > Og mixing is generated by 
two-loop mixing of O2 into O9 as well as by one-loop mixing and one-loop renormalization of 
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-^i,u-i4quark(^7)tree + -^i,J-i4quark(^9) 



the Qs factor in the definition of the operator Og. We denote the corresponding contribution 
to the counterterm form factors by Pf^'^^g ^^'^ ^Tu^^- obtain 
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The total counterterms F^^ (i = 1,2; j = 7,9), which renormahze diagrams la)-le), 
are given by 

pCtQ) _ pCtO) , pCt(j) , pCt(j) /Q^X 

Exphcitly they read 
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The quantities F^^^ (i — 1,2; j — 7,9) compensate the divergent parts of the form factors 
associated with the virtual corrections to 0"2- They are given by 



2,«,div g;^g' 

(29) 

<div = - 7^ [71820 + 21420 iTT - 252 5 - 275^ - 4.^^] - ^(32 - 17 L.), 



46 
243e' 

As mentioned before, we will take diagram If) into account only in Section IV. The same 
holds for the counterterms associated with the b- and ci-quark wave function renormalization 
and, as stated earlier in this subsection, the 0{as) correction to the matrix element of SZigOg. 
The sum of these contributions is 

SZ^{Oi) i.ioop + ^ ~ [<^^v(<^9)tree + (09)l-loop] , = Z^{mb)Z^{md) - 1, 

and provides the counterterm that renormalizes diagram If). We use on-shell renormal- 
ization for the external b- and o?-quark. In this scheme the field strength renormalization 
constants are given by 

Z4m)^l-^U^y(- + - + 4). (30) 

^ 47r 3 Vm/ \e eiR / 

So far, we have discussed the counterterms which renormalize the 0{as) corrected matrix 
elements {d£~^£~\0^\b) {i — 1,2). The corresponding one-loop matrix elements [of 0{a^)] 
are renormalized by adding the counterterms 

ajg ,^ V 
4 TT e 

F. Renormalized form factors of O'^ and O2 

We now have all ingredients necessary to present the renormalized form factors associ- 
ated with the operators O" and O^- We stress again that only the contributions of the 
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l,ti, div 



diagrams la)-e) and the counterterms discussed in Subsection HIE are accounted for in the 
result below. Diagram If) and associated counterterms will be included in the discussion 
of the virtual corrections to Og. We decompose the renormalized matrix elements of Oj 

{i = 1,2) as 

{d£+£-\Onb) = (09)tree + i^S(07)tree] , (31) 

where the operators O7 and Og are defined in Eq. (8). The renormalized form factors read: 



833 208 . AOiTT / 2 58 , 2i7r Am ^ Stt' 



l.U rron OAO ^ O/IQ I 'TOO O /I Q * O'/ O'Z * 



729 243 243 V 729 243 27 27 243 

/ 1453 14 ^ 2i7r Ain , Mtt^X „ 
+ ^ L. \ H ] (32) 

V 3645 243 ' 27 27 ' 243 y ^ ^ 

/_4712 68 2 2 _2i7r 267r2\ .3 

^ lv~5103^ 243 ^+27 ^~~9~^ ^ ' 

(9) _ 1736 224 2864 272 256 . 

520 in QAiTT ^ 272 m ^ 200 tt^ 256 ,,, 

\ Ls H \ C(3) 

243 243 243 729 27 ^ 

/ 388 20 , 16 , Ain Sin , 87r^\ ^ 
+ \ L.-\ L„ H L. s (33) 

V 675 27 ' 1215 27 27 ' 243 / ^ ^ 

/1018057 4 4 2 4 Ain 8n'^\ 

^ V 1786050 ^ 243 ' ~ 27 * ^ 2835 "^"9 81 J ^ 

92876363 344 ^ 4 ,2 16 ^ 20i7r 16i7r , 164 vr^^ 



48223350 729 9 ' 76545 27 27 729 

/7) 1666 416 , 80i7r / 4 116 , Am 8m ^ 16 tt^ 

Fk\] = \ L„ H h \ Ls 

243 81 ^ 81 V 243 81 9 9 81 



+ 



/2906 28 ^ Am 8m ^ 28 tt^ 



+ (34) 

V1215 81 9 9 81 / ^ ' 



/9424 136 4 2 4i7r 527r2\ ^3 
+ ivl70T""8r^^"9^^ + ^ stJ ' ' 

(9) _ 380 304 3136 544 512 . 

608 iTT 128 ivr 544 ivr 400 tt^ 64 

+ ^1 8i~^^ + ~8]~^'^"^43~ + "9"^^^^ 

/776 40 , 32 ^ 8m 16i7r , levr^X ^ 

+ Ls La \ Ls H s (35) 

V225 9 ' 405 ^ 9 9 ' 81 y ^ ' 

f 1018057 8 8 2 8 . SiTT 16n^\ 

+ 297675 ~81^+9^~ 945 ''"^^^J ^ 

/ 92876363 688 , 8 ,2 32 ^ 40 32i7r , 3287r2\ .0 

+ \ Ls + - L: La Ls H s^, 

\ 8037225 243 3 ' 25515 9 9 243 / ' 
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with 

Lg — ln(s) and = In 

As has been mentioned before, we only include terms up to 0{s^) in the result. We checked, 
however, that the terms of order are numerically negligible. 




IV. VIRTUAL CORRECTIONS TO THE MATRIX ELEMENTS OF THE OPER- 
ATORS O7, Og, O9 AND Oio 

The virtual corrections to the matrix elements of Oj, Og, O9 and Oiq and their rcnor- 
malization are discussed in detail in Refs. [34, 50]. For completeness we list the results of 
the renormalized matrix elements. They may all be decomposed according to 



as 



where 



47r 



.(1) 



(0) 



.(1) 



a 



(0) 



47r 



a. 



^ 47r 



a. 



.(1) 



and C?S=C^ 



A. Renormalized matrix element of O7 



The renormalized corrections to the form factors F^^^ and F^'' are given by 



^(9) 



16 A 1 . 1 .0 1 



(36) 



„C7'i 32 ^ 32 „ ^ „ 128 „o 
The function /inf collects the infrared- and coUinear singular parts: 



inf 



2e 



- l + s + -s^ + -s'^ 



1 2e 



(37) 



^ln(r) + ^ln(r)-^ln2(r), (38) 



eiR 3 V 2 3 

where eiR and r = {my ml) regularize the infrared- and coUinear singularities, respectively. 
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B. Renormalized matrix element of the operator Os 



The renormalized corrections to the form factors of the matrix element of Os are 
^(9) _ 104 32 2 40 2V ^4212 32 2V2 



/ 193444 560 A 16 . . . 
^(^)_ 32,^8 2 44_8_^/4, 40\ ^/32, 316\2 



^ 9'' 27 99 V3 3/19 9' ^ ^ 



/2OO 658\ 8 
V 27 9 y 9 ^ ^ 

C. Renormalized matrix element of Og and Oio 

The renormalized matrix elements of Og and Oio, finally, are described by the form factors 

(41) 

(42) 

(43) 
(44) 

where /inf is defined in Eq. (38). 

The contribution of the renormalized diagrams If), which have been omitted so far, is 
properly included by modifying Cg^^ as follows: 

For s < iz {z = ml/ml) the loop function Ho{z) can be expanded in terms of s/(4^). We 
give the expansion of Ho{z) as well as the result for Hq{0): 
1 





16 

= T + 


20 . 


16 .2 


116 .3 , 






H 






-^10 










-^10 











Ho{z) 



2835 



-1260 + 2520 In ( ^-^j + 1008 +432 T^^j + 256 



(45) 



/ X 8 4 , 4i7r 8 ^ 

ff.(0) = ^-glnW + — + -V 
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V. CORRECTIONS TO THE DECAY WIDTH B Xd£+£- 



The decay width differential in s can be written as 
dr{b 



ds 




l + 2s) 



+4(1 + 2/s) 



487r3 ' 

^pBrems, B 



CeflF 
9 



+ 



cf 



+^?E^ + ^1_. (46) 

as as 

The last two terms in Eq. (46) correspond to certain finite bremsstrahlung contributions 
specified in Appendix E. Their result can also be found in this appendix. All other correc- 
tions have been absorbed into the effective Wilson coefficients Cf^, Cf^ and Cfg. We follow 
[32, 34, 50] and write the effective Wilson coefficients as 



G 



eff 



l + ^u;Js] 



{A, - I Tga h{z, s)-^ h{0, s) + h{z, s) 



+ Ugh{l,s) + Wgh{0,s)) 



+ 



Cf 



1 + ^uj.is) ) 



(47) 



^10 



1 + ^u;,{s) ) Aio, 



TT 



where we have provided the necessary modification to account for the CKM structure of 
b — > dt^l~ . The renormalized form factors -Fi and Pl^jf' and can be found in Section III F 
while the renormalized form factors F^f\ F^f and F^'^^ arc given in [34, 50]. The 
functions u!7{s) and ujq{s) encapsulate the interference between the tree-level and the one- 
loop matrix elements of O7 and Og^io and the corresponding bremsstrahlung corrections, 
which cancel the infrared- and coUinear divergences appearing in the virtual corrections. 
When calculating the decay width (46), we retain only terms linear in ag (and thus in 
LOT, LOg) in the expressions for |C|^p and l^ioP- Accordingly, we drop terms of 
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0{a^) in the interference term Re (^Cf^Cf^*^ , too, where by construction one has to make 
the replacements Ug ujjg and ujj ujjg in this term. The function ug has aheady been 
calculated in [32], where also the exact expression for h{s, z) can be found. For the functions 
uj-j and ujjg and more information on the cancellation of infrared- and coUinear divergences 
we refer to [34]. 

The auxiliary quantities ^7, Ag, Aio, Tga, Tgb, Ug and Wg are the following linear combi- 
nations of the Wilson coefficients Cidi): 



As 
Ag 



^C.(,)-lc3(.)-lc.(,)-fa(,)-fCe(.), 

Csii^) + CM - \ CM + 20 CM - ^ CM: 



Att 4 

^-Cg{^) + -Cs{l^) 

asW 3 



64 64 



H9 



i=3 



In 



^10 

Tga 
Tgb 
Ug 
Wg 



^CM + C2(l^), 

■ 6 C3{ii) + 60 CM, 

7 2 32 

: - - C3(/x) - - CM - 38 CM - Y 

1 2 32 

: - - CM - 3 CM - 8 C^5(/^) - y C6(/.). 



(48) 



In these definitions we also include some diagrams induced by 03,4,5,6 insertions, viz the 
0{a^) contributions, the diagrams of topology If) and those bremsstrahlung diagrams where 
the gluon is emitted from the b- or d-quark line (cf [35]). 

For completeness, we give in Table I numerical values for Ci, C2, A7, As, Ag, Aiq, Tga, Tg^, 
Ug and l^g at three different values of the renormalization scale We note that the recently 
calculated contributions [33] to the anomalous dimension matrix which correspond to the 
three-loop mixings of the four-quark operators into Og have been included by adopting the 
procedure described in Appendix C. As can be seen in Table 1, some of the entries have a 
very small amount of significant digits. In our numerical analysis presented in Section VI 
we work with a much higher accuracy. 
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II = 2.5 GeV 



H = 5 GeV 



fj, = 10 GeV 



0.267 0.215 0.180 





1 / 


(-0.696, 0.240) 


(-0.486, 0.206) 


(-0.326, 0.184) 


(Cf , 




(1.046, -0.276) 


(1.023, -0.017) 


(1.011, -0.010) 






(-0.360, 0.032) 


(-0.321, 0.018) 


(-0.287, 0.009) 


\ 8 ' 


8 / 


(-0.169, - 0.015) 


(-0.153, - 0.013) 


(-0.140, -0.012) 




^9 J 


(^4.z4i, — u.uyi j 


(^4.izo, U.UDDj 


[4:.L6i, u.iyuj 


(rp(0) 




(0.118, 0.292) 


(0.376, 0.258) 


(0.577, 0.235) 






(-0.003, - 0.013) 


(-0.001, - 0.007) 


(0.000, - 0.004) 






(0.045, 0.023) 


(0.033, 0.015) 


(0.022, 0.010) 






(0.044, 0.016) 


(0.032, 0.012) 


(0.022, 0.009) 


US, 




(-4.373, 0.135) 


(-4.373, 0.135) 


(-4.373, 0.135) 



TABLE I: Coefficients appearing in Eq. (48) for /i = 2.5 GeV, fi = 5 GeV and ^ = 10 GeV. For 
as{n) (in the MS scheme) we used the two-loop expression with 5 flavors and as{mz) = 0.119. The 
entries correspond to the pole top quark mass rrit = 174 GeV. The matching for the top and for 
the charm contribution was performed at a scale of 120 GeV and 80 GeV, respectively [32]. The 
superscript (0) refers to lowest order quantities, while the superscript (1) denotes the correction 
terms of order a^, i.e. X = + X^^^ with X = C, A, T, U, W. Note that the contributions 
calculated recently in Ref. [33] are included. These contributions only affect the entries for a'^^ . 

VI. PHENOMENOLOGICAL ANALYSIS 

As the main point of this paper is the calculation of the NNLL corrections to the process 
b — > Xai'^i~, we keep the phenomenological analysis rather short. In the following we 
investigate the impact of the NNLL corrections on three observables: the branching ratio, 
the CP asymmetry and the normalized forward-backward asymmetry. As our main point 
is to illustrate the differences between NLL and NNLL results, we do not include power 
corrections (and/or effects from resonances), postponing this to future studies. 

Since the decay width given in Eq. (46) suffers from a large uncertainty due to the factor 
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m^pQig, we follow common practice and introduce the ratio 



R, 



■quark 



is) 



ds ^ ds 



(49) 



in which the factor m^p^j^ drops out. Note that we define -Rquarkl-s) as a charge- conjugate 
average as this is likely to be the first quantity measured. The expression for the semileptonic 
decay width r(6 — > XcePg) is as follows: 



r(6 



c,pole 



,pole 



mt 



(50) 



where 5^(2;) = 1 — 8 2; + 8 2;^ — 2;'^ — 12 2;^ ln(2;) is the phase space factor, and 



K{z) 



2as{mb) f{z) 



(51) 



incorporates the next-to-leading QCD correction to the semileptonic decay. The function 
f{z) has been given analytically in Ref. [51]: 



f{z) = -(1 - z') 



25 239 



25 



+ z ]n{z) 20 + 90^ 



4 



■z' + 
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+z' ln\z) (36 + z') + {l-z')(^^-fz + f z'^ ln(l - .) 

-4 (1 + 30 + /) In(^) hi(l - 2) - (1 + 16 + z^) (6 Li(z) - tt^) 

1 - 



-32 ^3/2(1 + ^) 



TT^ - 4 Li( V^) + 4 Li(- V^) - 2 In(^) In 



(52) 



ne- 



,1 + A/i 

In the following analysis we write the CKM parameters appearing in 6 — X^i'^i^ as 
glecting terms of C(A^)) 

with p — p{l — A^/2) and fj — 77(1 — A^/2) [52]. For Vcb, appearing in the semileptonic decay 
width, we use Vcb = AX^. Numerically, we set A = 0.81, A = 0.22, p = 0.22 and fj = 0.35. 
For the other input parameters we use as{mz) = 0.119, m^°^^ — 174 GeV, CKem — 1/133, 
rub = 4.8 GeV, mjnib = 0.29, mw = 80.41 GeV, mz = 91.19 GeV, and sin2(%) = 0.231. 

In Fig. 2 we show the /i-dependence of -Rquarkl-s) for 0.05 < s < 0.25. The solid lines 
correspond to the NNLL results, whereas the dashed lines represent the NLL results. We 
see that, going from NLL to NNLL precision, -Rquark(s) is decreased throughout the entire 
region by about 20 — 30%. Although the absolute uncertainty due to the //-dependence 
decreases as well, the relative error remains roughly the same. 
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FIG. 2: i?quark(s) as defined in Eq. (49). The solid lines show the NNLL result for /x = 
2.5, 5.0, 10.0 GeV, whereas the dashed lines show the corresponding result in the NLL approx- 
imation. At s = 0.25 the highest (lowest) curve correspond to /x = 10 GeV (/x =2.5 GeV) both in 
the NLL and NNLL case. 

As mentioned already in the introduction, the region 0.05 < s < 0.25 is free of resonances, 
as it lies below the J/^ threshold and above the p and uo resonances. The contribution of 
this region to the decay width (normalized by r(6 ^ XceUe)) is therefore well approximated 
by integrating i?quark('S) over this interval. At NNLL precision, we get 

0.25 

i?qua.k = j dsR^^^i,{s) = (4.75 ±0.25) x lO"^. (53) 

0.05 

The error is obtained by varying the scale between 2.5 GeV and 10 GeV. The correspond- 
ing result in NLL precision is i?quark = (6.29 ± 0.21) x 10"'^. The rcnormalization scale 
dependence therefore increases from ~ ±3.4% to ~ ±5.3%. The reason for this increace 
can be understood from Fig. 2: While for 0.13 < s < 0.25 the /i dependence of Rqaaxk{s) 
at NNLL and NLL precision is similar, the fi dependence almost cancels in the NLL case 
when integrating s between 0.05 and 0.13 due to the crossing of the dashed hues in this 
interval. This cancellation does not happen in the NNLL case, leading to a shghtly larger 
//—dependence of i?quark at NNLL. 
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-5 



0.05 




0.25 



s 



FIG. 3: CP asymmetry: The solid lines show the NNLL result for fx = 2.5, 5.0, 10.0 GeV, whereas 
the dashed lines show the corresponding result in the NLL approximation. At s = 0.25 the highest 
(lowest) curve correspond to /x = 10 GeV (/x =2.5 GeV) both in the NLL and NNLL case. 

As pointed out already, in the process b — > Xfii'^i" the contribution of the w-quark 
running in the fermion loop is, in contrast to b ^ Xsi~^i~, not Cabibbo-suppressed. As 
a consequence, CP violation effects are much larger in 6 — > Xdi^£~. The CP asymmetry 
acp(s) is defined as 



In Fig. 3 we show acp(s) for 0.05 < s < 0.25. The solid and dashed Unes correspond 
to the NNLL and NLL results, respectively. We find several differences between the two 
results: The sohd lines are much closer together. Also they cross each other at s 0.11. 
Furthermore, the NLL result clearly shows a positive CP asymmetry throughout the entire 
s region considered, while the NNLL lines indicate that acp(s) can be both positive and 
negative, depending on the value of s. Because of that, it does not make much sense to 
quantify the relative error due to the /^-dependence. The plot, however, clearly shows that 
the absolute uncertainty is much smaller in the NNLL result. For NLL results, see also 



dr{b^Xat+e-) _ dr{b^x^e+e-) 



acp(s) 



ds dS 



(54) 



dr{b-^Xde+e-) , dr{b-^x^e+e-) ' 



ds ds 



Ref. [53]. 
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We also give the averaged CP asymmetry ocp in the region 0.05 < s < 0.25, defined as 



0.25 



acp 



0.05 



ds 



ds 



0.25 



Jds I 



(55) 



0.05 



ds 



+ 



ds 



Varying between 2.5 GeV and 10 GeV one obtains the ranges 



1.4% < ag^L < 7.7%, ; 0.56% < ag^^^ < 0.93% . 



NNLL 



We now turn to the forward-backward asymmetry. As for i?quark(s), we introduce a 
CP-averaged version of the normahzed forward-backward asymmetry, defined as 



J d{cos9) sgn(cos^) 



d'^Tib-^Xdi+e-) , d'^r(b^X^£+£-) 



ds d{cos 9) 



+ 



ds d{cos 6) 



dr(b^Xai+i-) dTjb^x^e+e-) 

ds ds 



(56) 



where 9 is the angle between the three-momenta of the positively charged lepton and the 
6-quark in the rest frame of the lepton pair. The result of the integrals in the numerator 
of Eq. (56) for the case b Xgi^i'^ can be found in [38]. The corresponding result for 
b — > X(i£'^£~ is, up to different CKM-structures, the same. 



i<; 




S 



I I I I I I L_l LJ I I I I I L 



0.05 0.1 0.15 0.2 0.25 
s 



FIG. 4: CP-averaged normalized forward-backward asymmetry. The solid lines show the NNLL 
result for /x = 2.5, 5.0, 10.0 GeV, whereas the dashed lines represent the corresponding result in 
the NLL approximation. 
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In Fig. 4 we illustrate the //-dependence of Afb(s) in the region < s < 0.25. Again, 
the solid and dashed lines represent the NNLL and the NLL results, respectively. The 
reduction of the //-dependence going from NLL to NNLL precision is striking: one can 
clearly distinguish the three dashed lines, whereas the NNLL lines are on top of each other 
throughout the region. The position Sq at which the forward-backward asymmetries vanish, 
is essentially free of uncertainties due to the variation of /t: we find Sq = 0.158 ± 0.001. 
To NLL precision we get s^^^ = 0.145 ± 0.020. 

As a last illustration, we show in Fig. 5 the dependence of Rquaxkis) on the matching scales. 
In all the previous plots we used a matching scale of 120 GeV for the top contribution and 
a matching scale of 80 GcV for the charm contribution. In Fig. 5 the solid line corresponds 
to this scheme, while the dashed line is obtained by matching both contributions at a scale 
of 80 GeV. The difference between the two schemes is between 2% and 4%. 

0.5 
0.4 

in 

I 

O 

^ 0.3 
0.2 

3 

0.1 


0.05 0.1 0.15 0.2 0.25 
s 

FIG. 5: i?quark(s) for /i = 5 GeV. The solid line corresponds to matching top and charm contri- 
butions at 120 GeV and 80 GeV, respectively. The dashed curve is obtained by matching both 
contributions at a scale of 80 GeV. 

VII. SUMMARY 

In this paper we presented the calculation of virtual and bremsstrahlung corrections of 
0{as) to the inclusive semileptonic decay b — > X^i'^i'. Genuinely new calculations were 
necessary to attain the virtual contributions of the operators and O2 ■ Some of the dia- 
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grams (in particular diagrams Id) turned out to be more involved than the corresponding 
diagrams for the c-quark contributions. We used dimension-shifting and integration-by-parts 
techniques to calculate them. The main result of this paper, namely the n-quark contribu- 
tions to the renormalized form factors -Fi^, -^in^ ^^'^ ^2^1^ is given in Section IIIF. 

We shortly discussed the numerical impact of our results on various observables in the 
region 0.05 < s < 0.25, which is known to be free of resonances. As an example, we found 
the improvement on the forward-backward asymmetry Afb(s) defined in Eq. (56) to be 
striking: the NNLL result is almost free of uncertainties due to the //-dependence. 
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APPENDIX A: CALCULATION TECHNIQUES 

1. Reducing tensor integrals with dimension-shifting techniques 

We follow Ref. [42] and derive a method that allows to express tensor integrals in D 
dimensions in terms of scalar integrals of higher dimensions. 

An arbitrary L loop tensor integral with internal and E external lines can be written 
as a linear combination of integrals of the form (suppressing Lorentz indices of G^^^) 

G'^iM. {<]) = / (n|^) n %.,m' . (ad 



j=i 1=1 



where 



^ L E 

n=l m=l 



ki and qj denote the loop and external momenta, respectively. The matrices of incidences 
of the diagram, cu and rj, have matrix elements cOijjTjij e {—1,0,1}. The quantities {s^} 
and {m^} denote a set of scalar invariants formed from the external momenta Qj and a set 
of squared masses of the internal particles, respectively. Generically, the exponents Ui are 
equal to 1. However, often two or more internal lines are equipped with the same propagator. 
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This may be taken into account by reducing N to N^^ < N, thus increasing some of the 
exponents i/i. 

Applying the integral representations 



1 



da a" ^ exp i a (k'^ — -\-ie) 



(A2) 



and 



d 



exp 



i{aj kj) 



(A3) 



ai=0 



allows us to easily perform the integration over the loop momenta by using the D dimensional 

Gaussian integration formula 



k exp 



{Ae + 2(pk)) 



D 


r ■ 21 




A _ 



We find the following parametric representation: 



4: in 



j=i I 



d 



oo OO u—1 

daj a/ 
X / ... / ^ exp 



1 9{aj)^j, 

iQ{{s,},a) 



N 



— i ^ ar {ml —ie) 



r=l 



(A4) 



i 

The quantities Si are scalar invariants involving the external momenta qi and the auxiliary 
momenta Oj. D(a) arises from the integral representations of the propagators: let k be the 
L-dimensional vector that consists of all four-momentum loop vectors. The product of all 



Pt^ can then be written as 

Kj ,mj 



N / N 



Yl dtti f{a) exp i (k'^Bk + {bk) + cj 



with /cj-independent quantities f{a), B, & and c. D{a) denotes the determinant of the Lx L 
matrix B. 

The differentiation of G^^^ in Eq. (A4) with respect to aj generates products of exter- 
nal momenta, metric tensors gf^j^ and polynomials R{a) and provides an additional factor 
D{a)~^. Because of 



R{a) exp 



N 

E 

r=l 



R{id) exp 



TV 

E 

r=l 



, with dj 



d 
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we may replace the polynomials R{a) with R{id). The additional factor of 1/D{a) can be 
absorbed by a redefinition of D, i.e. by shifting D to D + 2 and multiplying with a factor 
(Ain)^. The crucial point is that all factors generated by differentiation with respect to aj 

may be written as operators which do not depend on the integral representations we have 
introduced in Eqs. (A2), (A3). Therefore, it is possible to write tensor integrals in momentum 
space in terms of scalar ones without direct appeal to the parametric representation (A4): 



N 



n 



N 



where the tensor operator T (suppressing its Lorentz indices) is given by 



(A5) 



=exp[-ig({si},Q;) (4i7r)^</+] 



X 



N rij 

nn 



d 



exp 



iQ{{si},a) {4i7rfd+ 



(A6) 



ai=0 



aj=iOj 

The operator cZ+ shifts the space-time dimension of the integral by two units: 

({.-.}, K}) = ({.-.}, 

Notice that throughout the derivation of the tensor operator T the masses mj must be kept 
as different parameters. They are set to their original values only in the very end. 



2. Integration by parts 

According to general rules of D dimensional integration, integrals of the form 

_d K 

vanish. There may exist suitable linear combinations 

■I *■ n,'L,(*i-™j+i.)''' 

that lead to recurrence relations connecting the original integral to simpler ones. The task 
of finding such recurrence relations, however, is in general a nontrivial one. A criterion for 
irreducibility of multi-loop Feynman integrals is presented in [43]. In [42], the method of 
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partial integration is combined directly with the technique of reducing tensor integrals by 
means of shifting the space-time dimension. 
The integral 

1/1^21^31^41^5 I LI b LI I 1^11/21/31/41/5 

"^""^ "^""^ [P] [r^] [{I + r)2] [(/ + qy] ^ [(r + p)' - ml] ^^^^ 

enters the calculation of diagrams Id). At the same time it is a very good example to 
illustrate the integration by parts method. The operators 1^, 2^,... are defined through 

i±r - T 

^ ^1/11/21/31/41/5 — -'t'lltl l/2t'3f'4t'6 ) ■ ■ ■ ■ 

The present case is especially simple because we only need to calculate one derivative. Using 
the shorthand notation ^1/11/21^31/41/5 = /{j,.} we get (for i/j > OVi): 

d 



r^I{ui} = -D-2z/2r^2+ -2z/3r(/ + r)3+ -2z/5r(r + p)5+ 



Scalar products of the form (a h) we write as [a? + 6^ — (a — 6)^] /2 and find 

^r^'hu,} = [d-2u^-u^-u^- i/3(2- - 1-) 3+ - 1/52-5^^ ^ 

At this stage we might also reduce some of the scalar products by shifting the dimension. 
The corresponding procedure is presented e.g. in [42]. In the present case, however, the 
pure integration by parts approach suffices. The identity 

d^'r J- = 



yields directly the desired recurrence relation for the integral 



U1U21/31/41/5 ^ _ 2 — 1/3 — 1/5 1^11^21/31^41/5' y-^'^j 



Subsequent application of this relation allows to express any integral F^^^{i'i} with indices 
Ui e N"*" as a sum over integrals F'^^^{i'i} with at least ui — or 1/2 — 0. 
The general procedure is the following: 

• One expresses suitable scalar products in the numerator of a given Feynman integrand 
in terms of inverse propagators Pj^^ and cancels them down. It is important to notice 
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that it is not always the best strategy to try to cancel down as many scalar products 
as possible. The resulting set of integrals to calculate highly depends on which scalar 
products one cancels down. The best way is to try a couple of different cancellation 
schemes and compare the resulting integrals. 

• One writes the integral as a sum over tensor integrals of the form (Al) with products 
of /cf . For each of those integrals the tensor operator T is determined in order to 
reduce the problem to scalar integrals with shifted space-time dimension. 

• One applies appropriate recurrence relations to reduce the number of propagators in 
the integrals, hoping to be able to solve the remaining integrals. 



APPENDIX B: CALCULATION OF THE DIAGRAMS Id) 

The contribution of the sum of diagrams Id) is given by a combination of integrals of the 
form 

[P] [t^] {{I + Tf\ {{I + qY\ {{t + V? - m?] ■ ^ ^ 

In this section we show how to solve these integrals with the methods presented in Appen- 
dices Al and A 2. The function £>(«), which is independent of ni and rir-, is not needed in 
order to find the tensor operators T. Nevertheless, we give it as an illustration: 

D{a) = {ai + as + 04) {a2 + + ^5) - aj. 

The function Q(^{si}, a) , however, must be calculated for each type of tensor integral. As 
an example we give Q{{si}, a} for ni = 0, = 1: 

Q{{si}, a) = — (cti + a3 + 04) (aip) — q;3Q;4 (aiq) — - (cti + as + "4) «i- 



The corresponding tensor operator T reads: 
TP'{q,p,d,d+) = 167r^d+ 



The action of an operator di on the integral F^^y is 
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The next step is to repeatedly apply the recurrence relation (A8) on the integrals Fv^XviVif^ 
until vi or becomes zero. The problem is then reduced to the calculation of the two types 
of integrals 

J, ,. 
In the present calculation D may take the values 



D = 4 - 2e, 6 - 2e, 8 - 2e or 10 - 2e. (B4) 

It is important to note here that the denominator in Eq. (AS) can become proportional to 
e for certain values of 1/2, J^s and 1/5. Thus, some of the integrals in (B3) need to be 
calculated up to C(e^). 

The first type of integrals (-^0,^131/41/5) can easily be solved individually by using a single 
Mellin-Barnes approach. This method naturally results in an expansion in s. Furthermore, 
the occasionally needed 0{e^) terms are easily obtained since the expansion in e is done only 
in the very end. We now turn to the much more complicated calculation of the second set 
of integrals. Instead of calculating every single occurring integral individually, we derive a 
general formula for -Fj^tsi^^i^s where we are left with a three-dimensional Feynman parameter 
integral: 

1 

K,o.3.,.,-{ r(i.i)r(i/3)r(i/4)r(i/5)7 ^ 



X(l - ^x)^/2-'^3-1^.3-l (1 _ ^Y3+.,-l-D/2yU,-l _ ^)D/2-l-.i 

A = (1 — x) (1 — uy) — sxyu — iS. 

We now replace all occurrences of -Fj^tsi^^^g according to Eq. (B5) and are left with a three- 
dimensional integral over a rather lengthy integrand. This integrand can be split up into 
three different parts: 

• A part with no additional divergences arising from the integrations. 

• A part with problematic ^-integration. 

• A part with problematic li-integration. 
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In the first part, the regulator e is not needed at all and may be set equal to zero at the 
very beginning. The occurring integrals can then either be performed directly or with the 
use of a single Mellin-Barnes representation. The second part boils down to two different 
integrals, which can both be computed with subtraction methods. The last part is clearly 
the most difficult one. It can be reduced to three integrals which we calculate using a 
double Mellin-Barnes representation. Since this double Mellin-Barnes is very different from 
the one presented in Subsection 3.1.4 in [34], we give, as an example, the needed procedure 
to calculate one of the three integrals. Specifically, we have to deal with the integrals 
111 



[1 — m)^+^ ((1 — x){l — uy) — sxuy — i5)' 



/, = idx jdy jdu ^^ . ....^^. ^ m 





where j can take the values 0, 1 and 2. We focus on the case where j = 0. We introduce a 
first Mellin-Barnes integral in the complex i-plane with the identifications (for notation see 
e.g. [34]): 

^ (l-x)(l-uy), sxuy + i5, X^l + 2e, 



and get 



111 

g-i7r(l+2e) 



g-i7r(l+2ej f r f f 



U^+'^y^{l — yY ^T* 
;i - Uy+'{1 - x)l+^+*(l - M2/)l+2^+* 



The path 7 lies in the left half-plane and can be chosen arbitrarily close to the imaginary 
i-axis. We introduce a second MeUin-Barnes representation in the complex f-plane for the 
last factor in the denominator of Eq. (B7). For this, we rewrite \ — uy d& 1 — u ^ u{\ — y) 
and make the following identifications: 



yielding 



K^^u{\-y), ^ -{\-u), A = l + 2e + i, 



111 

g-i7r(l+2e) 



h = (2z7r)^ r(l + 2e) ' ' ' ' \ duV{-t)V[-i!) 



7' 7 



iM\ — nA^-Se-t-t't 
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The path 7' hes to the left of the imaginary i'-axis and can again be chosen arbitrarily close 
to that axis. The parameter integrals can now be performed and give products of Euler 
Beta-functions. We work out the remaining integrals over t and t' applying the residue 
theorem. For this, we close the t-integral in the right half-plane and focus on the enclosed 
poles. There are two different sequences of poles, namely poles that depend on t' (coupled 
poles) and poles that do not (uncoupled poles). The latter poles lie at the following positions: 

• t^0,l,2, S, ... , 

• t — — e, 1 — e, 2 — e, . . . , S — €, . . . , 

Note here that Jq exists only for negative values of e. The pole located at t = — e therefore 
lies in the right half-plane and needs to be taken into account. Since we are interested in an 
expansion in s, we can truncate the two pole sequences at a suitable S. After calculating 
the necessary residues, we close the ^'-integral in the right half-plane as well and are arrive 
at pole sequences situated at the following positions: 

. f = 0, 1, 2, 

• t' = - e, 1 - e, 2 - e, . . . , 

• = 2 - - 3e, 3 - - 3e, 4 - - 3e, . . . ior t ^ N, N eN, 

t' = 2 - N - 2e, 3 - N - 2e, 4 - N - 2e, . . . for t = N - e, N eN. 

For A" > 3, some of the poles above he in the left f-half-plane and must be omitted. Unlike 
the procedure given in Subsection 3.1.4 of [34], we need to sum up the residues of all poles 
in the enclosed area. 

Calculating the contributions of the coupled poles in t, which lie at t = 2+n— 3e— t', n G N, 
yields an expression that is proportional to 52+n-3e-t' 'j^q problems now arise if one closes 
the integration path of the ^'-integral in the right half-plane: due to the —t' in the exponent 
of s, one gets an expansion in inverse powers of s, forcing one to calculate the residues of 
all enclosed poles. The second problem is even worse: for any given value of n, there always 
exists an infinite pole series in t' which contributes to the desired result. Thus, one also 
has to consider the infinite pole series in t. In order to avoid these problems, we close the 
integration path in the left half-plane of t'. The poles are then located at 

• t' — —1 — e, —2 — e, —3 — e, ... , 
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• t' ^-1- 2e, -2 - 2e, -3 - 2e, . . . , 

• f = -1 - 3e, -2 - 3e, -3 - 3e, . . . . 

After calculating the necessary residues we obtain the result for /q- The results for Ii and 
I2 are calculated in an analogous way. 

APPENDIX C: SOLUTION OF THE RENORMALIZATION GROUP EQUA- 
TION FOR THE WILSON COEFFICIENTS 

The Wilson coefficients satisfy the renormalization group equation 

= 7^(«s) C{f,) , (CI) 

where 7(q;s) is the anomalous dimension matrix. This matrix can be written as a Taylor 
series in ccg: 

The general solution of Eq. (CI) can be expressed with the evolution matrix C/(/x, /xq): 

[/(/io,/io) = 1. (C2) 
The aim in this section is to find a handy expression for U{ii, /Iq). 

The matrix 7^°) can be diagonalized. We introduce new quantities in the following way: 

(7(//) = veil,), 

^ i ^ 0,1,2,... , (C3) 

U{i^,l^o) = VU{fi,fio)V-\ 

The matrix V is chosen such that 7^°^ is diagonal. One can check that the new quantities 
satisfy equations similar to (CI) and (C2): 

^g(,)^7>.)g(.), 

Cifx) = [/(//, //o)C(//o), (C4) 
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We will now construct a solution to Eq. (C4). Once this solution is found, we can easily gain 
the solution of the initial problem for the non-diagonal 7^°^. The evolution matrix C/(/x, //q) 

— * 

satisfies the same equation as itself: 

^^C/(/x,/xo) = 7^(q;s)C^(/^,/^o)- (C5) 
We make the following ansatz for [/(//, //q): 

U{^^, /.o) = (1 + E (^) ^ ^) (/^' /^o) K, (C6) 
where C/^°^ (/x, //q) solves Eq. (C5) to leading logarithmic approximation and is given by 



D 



The vector 7*^°) collects the diagonal elements of 7*^'^^ The matrix K must be chosen such 
that the boundary condition given in Eq. (C4) is met. The quantities Pi, i — 0,1,2,... 
appear in the RGE for ccg: 



i+2 



Inserting the ansatz (C6) into Eq. (C5) and using the exphcit expression for C/(°)(/x, //q), the 
Ihs and the rhs of this equation can be written as 

j=i V ^ / 
j=i \ / 



The unknown matrices Jj can now be constructed order by order in ctg through the relations 
Lj = Rj. We give the explicit solutions to Ji and J2 since we need them to find the Wilson 
coefficients Ci{ii) to NNLL precision: 



7 — X -(o)j^i % (r'7\ 

'Jl,ij — 'Jijli r,ri2 ^r, ^fO) ^fOl ' y^^> 



Tr+(2A-|7f)^«.+ (7<"'^.)« 



J2,i3 - Sij% _ . ^(0) Z(0) (^8) 



36 



The result for Ji agrees with the one given in Section III of [45] . After we did the calculation 
for J2, we found out that the result already exists in the literature [54]. The two results 
agree as well. 

The matrix K is given through 



J,-ji)-rO{al). 



4. - V 4. 

With these informations at hand, we can present the evolution matrix for the initial problem 
given in Eqs. (CI) and (C2): 

[/(/.,/.o) = V \ l + ^^J,+ {^^^' j)iU^'\t,,t,o)KV-' + 0{al). (CIO) 
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APPENDIX D: ONE-LOOP MATRIX ELEMENTS OF THE FOUR-QUARK OP- 
ERATORS 



In order to fix the counterterms Fl^SSiq\xa.rk = 1^2) in Eq. (22), we need the one- loop 
matrix elements {d l'^ l~ \0 j\h) i.\oop of the four-quark operators O", O2, O4, and 0"2- 
Due to the 1/e factor in Eq. (22) they are needed up to C(e^). The explicit results read 

^2 + 3i7r-3Lj + 



4 4 

9e 27 



e 

81 



52 + 24 iTT - 21 TT^ - (24 + 36 m)Ls + 18 \ {09)tree , 



(d^+r|04|&)i-ioop = - ( — 



2e 



1 + ^(705 + 75^ + 5^) 
9 945 ^ ^ 



(O7) 



tree 



+ 



16 



27 e 8505 
4e 



-420 + 1260 in - 1260 + 252 s + 27 + 4 s^) 



+ (420 in + 910 - 630 in - 420 - 315 tt^ 

8505 



+ 315 - 126 s + s^) 



9/ tree 



37 



{d tr 1 0^2 1 b) i-ioop = I (c? i^r 1 1 b) i.ioop 



APPENDIX E: FINITE BREMSSTRAHLUNG CORRECTIONS 

In Section V those bremsstrahlung contributions were taken into account which generate 
infrared and collincar singularities. Combined with virtual contributions which also suffer 
from such singularities, a finite result was obtained. In this appendix we discuss the remain- 
ing finite bremsstrahlung corrections which are encoded in the last two terms of Eq. (46). 
Being finite, these terms can be directly calculated in d — A dimensions. 

The sum of the bremsstrahlung contributions from O7 — Og and Og — O9 interference 
terms and the Os — Os term can be written as 

^pBrems, A ^pBrems ^pBrems ^pBrems 

78 _|_ 89 _|_ 88 

ds ds ds ds 

J li^j 487r3 X (2Re[c78T78 + C89r89]+C88T88), (El) 



where 

2 



C78 — ■ <-^7 <-^8 ' "^89 — <-^F ■ <-^8 '-^9 ' "^88 — 



(0,eflf) 



(E2) 



For the quantities r78, t^^ and r88 we refer to [35]. 

The remaining bremsstrahlung contributions all involve the diagrams with an or 
2 insertion where the gluon is emitted from the u- or c-quark loop, respectively. The 

corresponding bremsstrahlung matrix elements depend on the functions Ai^^''^\ In 

d — A: dimensions we find 

xy(l -yf 



a4? = 8 (5 r) j\x dy ^^^^^ , Aig = 8 (? r) j^x dy : 



4? = 8 (g r) dy ^^(^ , Aig = 8 (g r) j\x dy 
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where 



C(") = -2xy{l- ij){qr) -q\J{l-^J)- z 5 , 
C'('^) — ml—2xy {1 — y){qr) — y {1 — y) — iS . 



The analytical expressions for Ai^23 Aigy can be written in terms of functions Gi{t) 
Azg = -2 + 



w — s 



where 2; = ml/ml. Gkit) {k > —1) is defined through the integral 



1 

Gkit) = Jdxx'' ln[l-tx{l- x) -i6], G^{t) = ^Go{t). 


Explicitly, the functions G-i(i) and Go{t) read 

2 TT arctan (y^) - t ~ ^ arctan^ (-^Z^) , ^ < 4 

C-iW = < 

^_2^^ ln(^±^) - ^ + 2 ln2(^*±|0) , ^ > 4 



7ry^-2-2y^arctan(^y^)), t<4 



G'o(i) = < 



^_^7ry^-2 + 2y^ln(^%^) , t >4 



The quantities Ai^"'* we obtain from Ai^'^'' in the limit 2; — > 0: 

A^(^) = -2 + -^[lnH-ln(5)], 



Ac) 



w — s 



A4? = -2[lnH -ln(s)] 



Following [35] , we write 



ds 



47r/ V47r 



487r3 



X 



1 

y (iu; I (cii + C12 + C22) T22 + 2 Re [ (ci7 + C27) T27 + (ci8 + 023) T28 + (cig + C29) T29] I ■ 
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Expressed in terms of the quantities Ai|f and Ai|f , defined by 



the quantities r^j introduced in Eq. (E7) read 
8 {w-s){l-wf 



T22 



27 



X i[3w^ + 2s2(2 + w)-sw(5-2w)] |Ai^ff + 



(E8) 
(E9) 



2s^(2+ + sw(l + 2w)] |Ai|^|V4s (1 --«;)- s (2 + w)] • Re [ Ai^f Ai^f ] I , 



(ElO) 



T27 = ^— X 
6sw 



{1 — w) {As^ — sw -\- w^) + s w (4 + s — w) ln( 



w 



At 



eff 
23 



4s^(l-w) + sw(4 + s-w) ln(w) Ai^B, (Ell) 



7"28 



9 siy (ly — s) 



X 



(w-s)2(2s- w)(l -w) 



A.^23 - 



2s(w - s)2(l - w] 



+ sw 



(1 + 2 s - 2 w) Ai^l - 2 (1 + s - w)M 



In 



(1 + s — w){w'^ -\- s{l — w)) 



(E12) 



4 1 



2 s[l — w){s + w) + 4 sw In(u') 

2 s(l — + + ^(3 s -\-w) \a.{w) 



Ai 



(E13) 



The coefficients cij include the dependence on the Wilson coefficients and the color factors. 



cii = C, 



Tl 



(0) 



C12 = C't2 ■ 2 Re 



(-y(O)^(O)* 



C22 — 



(0) 



^ ^(0)p;(0,eff)* 
Ci7 — Ot-j • O7 , 

^ _r r^(o)p;(o,cfr)* 
^ ^(o)p;(o,efr)* 

C19 — *-^T2 ■ *-^9 ) 



^ ^(0)P;(0,eff)* 
C27 — Up ■ ^2 ^7 1 

C28 — ■ L^2 '-^S ! 
C29 — Of • O2 O9 



(E14) 
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The Wilson coefficients C^jl g are given in Eq. (47) and numerical values for the coefficients 
Cf^ can be found in Table I. The color factors Cp, Cr-^ and C^-j are presented in Section II. 
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